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Abstract. We present a review of the normal form theory for weakly dis- 
persive nonlinear wave equations where the leading order phenomena can be 
described by the KdV equation. This is an infinite dimensional extension of 
the well-known Poincare-Dulac normal form theory for ordinary differential 
equations. We also provide a detailed analysis of the interaction problem of 
solitary waves as an important application of the normal form theory. Sev- 
eral explicit examples are discussed based on the normal form theory, and the 
results are compared with their numerical simulations. Those examples in- 
clude the ion acoustic wave equation, the Boussinesq equation as a model of 
the shallow water waves, the regularized long wave equation and the Hirota 
bilinear equation having a 7th order linear dispersion. 
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1. Introduction 

In this chapter, we review the normal form theory developed in |l4|, |l5| for 
weakly nonlinear and weakly dispersive wave equations where the leading order 
equation is given by the KdV equation in an asymptotic perturbation sense. The 
chapter is based on the report |l6| and the master thesis of the first author at Osaka 
University. 

A series of lectures was carried out by the second author in the Euro Summer 
School 2001 held at the Isaac Newton Institute, Cambridge for August 12-25, 2002. 
The lecture started with a brief summary of the Poincare-Dulac normal form theory 
for a system of ordinary differential equations ||. The main point in the lecture 
is to present the normal form theory for near integrable system where the leading 
order system is given by a nolinear wave equation. This may be considered as 
an infinite dimensional extension of the Poincare-Dulac normal form theory. The 
basic technique of the normal form theory based on the Lie transformation can be 
formally extended to the infinite dimensional case. Unfortunately the convergence 
theorem of the normal form series has no extension to the present theory. However 
one should emphasize that the leading order in the present theory is given by a 
nonlinear equation, and it is not clear how one define a resonant surface for the 
leading order equation. This may provide a good problem for a future project. 

Let us briefly summarize a background of the normal form theory for near inte- 
grable systems of nonlinear dispersive equations. It is well-known that for a wide 
class of nonlinear dispersive wave equations, the leading order nonlinear equation 
in an asymptotic expansion turns out to be given by an integrable system, such as 
the Korteweg-de Vries (KdV) equation in weak dispersion limit and the nonlinear 
Schrodinger equation in strong dispersion limit (see for example [^8) ) . This implies 
that most of the nonlinear dispersive wave equations are integrable at the nontrivial 
leading orders in an asymptotic sense. Then a natural question is to ask how the 
higher order corrections affect to the integrability of the leading order equations. 
In |l9) , the effect of the higher order corrections on one-soliton solution of the KdV 
equation was studied, and it was shown that the velocity of soliton is shifted by 
the secular terms in the higher order terms. Those secular terms or resonant terms 
are given by the symmetries of the KdV equation. The nonsecular terms then con- 
tribute to modify the shape of soliton. However the multi-soliton interactions were 
not studied in that paper. In p^ j, the normal form for weakly dispersive equations 
was first introduced up to the second oreder corrections, and it was found that the 
integrable approximation can be extended beyond the KdV approximation but not 
to the second order. This obstacle to the asymptotic integrability plays no rule for 
one-soliton solution, but provides a crucial effect for two soliton interactions. This 
was found in jl5], |l6) . The obstacles are defined as the nonexistence of the integrals 
of perturbed equation in the form of the power series in a small parameter. This 
was also recognized as the nonexistence of approximate symmetries of the perturbed 
equation [pl| . (The normal form for strongly dispersive wave equations has been 
also studied in 0, (jJJ.) The n in (ll| , the effect of the obstacle on the interaction of 
two soliary waves was studied for the regularized long wave equation (although the 
method can apply to other equations of weakly dispersive system). An inelasticity 
due to the obstacle was found, and it leads to the shifts of the soliton parameters 
and a generation of a new soliton as well as radiations through the interaction. 
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In this lecture, we present a comprehensive study of the normal form theory for 
weakly dispersive wave equations: We start in Section 2 to define the perturbed 
KdV equation as an asymptotic expansion of a weakly dispersive wave equation 
whose leading order equation is given by the KdV equation. We give a recursion 
formula to generate the higher order corrections which may be obtained by an 
asymptotic perturbation method (see for example p8|). The set of those higher 
order terms forms an extended space of differential polynomials which includes some 
nonlocal terms. The space is denoted by P Q dd, where "odd" implies the weight of 
the polynomial. 

In Section 3, the conserved quantities or integrals of the KdV equation are re- 
viewed, and we discussed approximate integrals of the perturbed KdV equation. 
We then obtain the conditions for the existence of approximate integrals in each 



order (Proposition 3.1). We also discuss a connection of the conserved quantities 
and the iV-soliton solutions of the KdV equation. 

In Section 4, we review the symmetries of the KdV equation, and discuss the 
approximate symmetries for the perturbed KdV equation. Here we also define the 
space of 'Pcven, which together with Podd provides the appropriate spaces for the 
normal form transform defined in the next section. 

In Section 5, we describe the normal form theory. The normal form transforma- 
tion is then obtained by a linear equation of an adjoint map defined as 

ad^(0) : Vcvcn ► Podd, 

where is the KdV vector field. The explicit form of the normal form is given 
for the perturbed KdV equation which contains the first three lowest weight ap- 
proximate conserved quantities (Proposition 5.1). The normal form then admits 
one-soliton solution of the KdV equation, which confirms the result in We 
also discuss the Gardner-Miura transformation which is an invertible version of the 
Miura transformation, and show that the inverse Gardner-Miura transformation is 
nothing but the normal form transformation after removing the symmetries of the 
KdV equation (Proposition 5.2). 

In Section 6, we consider the interaction problem of two solitary waves, and 
provides the explicit formulae for the shifts of the soliton parameters (Proposition 
6.1). we also give the formulae of the radiation energy and additional phase shifts 
which are used to compare with the numerical simulations for some examples in 
the next section. 

Finally in Section 7, we present explicit examples including ion acoustic wave 
equation, the Boussinesq equation as a model of shallow water waves and the regu- 
larized long wave equation. We show the good agreements with the results obtained 
by the normal form theory. We also consider a 7th order Hirota bilinear equation 
which admits an exact two soliton solution but is known to be nonintegrable. We 
look for an obstacle to the asymptotic integrability, and find that the obstacles 
appear at the fourth order. This implies that the forth order obstacles play no rule 
for two soliton solution, just like all obstacles has no rule for one-soliton solution 
for the system with the first three approximate integrals. 



2. Perturbed KdV equation 



Under the assumption of weak nonlinearity and weak dispersion, the wave prop- 
agation in a one-dimensional medium can be described by the KdV equation in 



4 



HIRAOKA AND KODAMA 



the leading order of an asymptotic expansion. Using an appropriate asymptotic 
perturbation method (e.g. see 28 ) one can show that the higher order correction 
to the KdV equation has the following expansion form with a small parameter e 
with < e < 1, 

(2.1) u t + K(u;e) = 0(e N+1 ) 1 

with K(u; e) = K<®(u) + eK^(u) + e 2 K^(u) + ■■■ + e N K^(u), 

where K^°\u) gives the KdV flow and the higher order corrections K^(u) are 
generated by a recursion formula starting from n = — 1, 

M(n) 

KM(u) = £ a^X^(u) 
( 2 - 2 ) _ >) 



af ] u [2n+3)x + £ (X^D-'X^) (u). 



1 +r, 2 = n-2 

i<»<m (m) 

l<J'<M(n 2 ) 

Here a J -™\ cj^ are the real constants determined by the original physical problem, 
u n x = d n u/dx n , and D _1 indicates an integral over x, Z? _1 (-) := f_ dx'{-). Each 

.Xj (u) is a monomial in the polynomial if^(u), and M(n) is the total number 
of independent monomials of the order n. The first few terms of K^ n \u) are then 
given by 

A''- 1 ' = a£ -1) u x> M(-l) = l, 

A-(°) = afW + <4°W, M(0)=2, 
A (1) = a^tiBa + a^u 3a; u + a^u 2x u x + a^u x u 2 , M (1) = 4, 

W21 ( 2 ) , (2) (2) (2) 2 , ( 2 ) 

Aa ' = a} 'u 7x + a 2 'u 5x u+a 3 'u ix u x + a\ 'u 3x u +%' u 3x u 2x 

+a^u 2x u x u+ a^Uxift+ajjul, M(2) = 8, 

ts(3) (3) , (3) i (3) . (3) 

K y ' = a\ u 9x + a 2 u 7x u + a 3 'u 6x u x + a\'u 5x u 2x 

, (3) 2 , (3) , (3) . (3) 

+ag 'u 5x u + a 6 'u4 X u 3x + a) 'u 4x u x u + u 3x u 2x u 

(3) 2 , (3) 3 , (3) 2 (3) 2 

+ag u 3x u x + a\ 'u 3x u° + a\{u 2x u x + a\ 2 'u 2x u x u* 

+a{ 3 3 ] u x ii + a ( ^u x u 4 + UxD' 1 (ul) , M(3) = 15. 

We normalize the KdV flow so that K^°'(u) is given by the standard form which 
we denote by Ag (u), i.e. 

Each polynomial A"(") has the scaling property: Assign the weight 2 to u(x, •), and 
1 to d/dx. Then if u(x. ■) = S 2 v(Sx, •) = (5 2 w(£, •), we have 

K^(u(x, •)) - S 2n+5 K (»>(„(£, •)) = ( 5 2n+5 (a( r ^ ( 2„ +5 )« +•••)■ 

Thus each polynomial K^ n \u) has the homogeneous weight "2n + 5". We denote 
by PoddM the set of all the odd weight polynomials generated by the formula (2.2). 
Then we have 

oo 

PoddM = ^ddM, 

n=-l 
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where pIjhM * s the finite dimensional subspace of the polynomials with the homo- 
geneous weight In + 5, and the dimension of the space is given by dim PiljdM = 
M(n), " 

P™[u] = Span R { *<">(«) : 1 < i < M(n) } . 

As shown in the examples, the subspaces PoddM up to n = 2 are given by the 
differential polynomial of u and its derivatives, and we denote 

^dM = fWsM ©r Q2„+ 5 M, 

where Pfefu] is the space of homogeneous differential polynomials of weight k, 

V k [u] = Span R J u l °u 1 ^ ■ ■ -ufa : ^(j + 2)1 s = k, I, e Z> 
and Qk [u] consists of the polynomials of weight k containing the integral operator 

Q, n+5 [u] C Span R [x^D^X^ • ■ ■ D^X™ : " = n I " 

The space 'P/Ju] can be also extended for k =even, and we will later define Qfc[u] 
for A: =even. 



Remark 2.1. In order to verify the expansion (2.1), one may need to impose the 
following conditions for the initial data on u(x,t = 0), 

a) |-u(a;,0)| < Cexp(— e 1 / 2 ^!), as \x\ — > oo, 

OO p 

b) \\u(x,0)\\ 2 HCO(K) / \u nx (x,0)\ 2 dx < oo. 

The regorous justification of the KdV equation from a physical model such as the 
shallow water waves has been discussed in [|[ [To| . 

3. Conserved quantities and N soliton solutions 

The integrability of the KdV equation implies the existence of an infinite number 
of conserved quantities or integrals. Here we briefly summarize those quantities, 



and discuss the approximate integrals for the perturbed equation (2.1). The ap- 
proximate integrals will play a fundamental role for the nomal form theory discussed 
in Section and provide an analyical tool to study the nonitegrable effect on the 
solution of the perturbed KdV equation. 

3.1. Conserved quantities. Let us first recall the definition of a conserved quan- 
tity for the evolution equation in the form, 

(3.1) Ut = f(u), with f(u) eP dd[u}. 



Definition 3.1. An integral of a differential polynomial p(u) S V[u] = ®kPk[u], 
I[u] = / p(u) dx, so that p{u) S V[u]/Im(D), 
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is a conserved quantity of ([O]) if 

p t e Im(-D), with D = d/dx. 
The polynomial p(u) is called a conserved density for (|3.l| ). 

We also define a vector field generated by /(it) and its action on the space V[u]. 

Definition 3.2. A vector field generated by f(u) is defined by 

d 



dm, 



i=0 

which acts on the space V[u] as a differential operator, 

V f : V[u] — ► V[u] 
9 1 — > Vf-g 

With this definition, the condition for p(u) being a conserved quantity can be 
expressed by 

p t = V f ■ p e Im(£»), 

The KdV equation with f(u) = K^(u) e V 5 [u] has an infinite sequence of con- 
served quantities, 



1 k 



' M = / pf ( u ) dx > with Pk ] e v 2k+2 [«] for fc = 0, 1, 2, • • • . 
which are generated by the bi-hamiltonian relation, 
(3.2) DVI^\(u) = eV/fV). with e = L> 3 + 2(Du + wD). 

The gradient (VI) (u) is defined by 

/* d 
I vCVI)(u) dx = lim —I[u + 6v], 
Jr S^odS 

which can be expressed as 

(vj)(u) = ^(-i^^-M for J M = / dx - 

i=Q ra ® 

The first few conserved densities pi are given by 

(0) 1 (0) 1 2 

Pa 9 U ' Pl 2 U ' 

pf ] = \{ul - 2u% = i« - 10-4 + 5« 4 ). 

Each density p^ (u) can be considered as 

pT eV 2k+2 [u}/Im(D). 

We then define 
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Definition 3.3. The set of differential polynomials for the conserved densities are 
given by 

F k [u]=V k {u]/Im(D), 

where ^[u] are defined by 



Pk[u] '■= Span R < u l °u l ^u' 



n 
3=0 



The following table shows the relation between the weight and the dimension of 
the space Tk [u] : 



Table 1 : The relation between the weight and the dimension of [u] 



Weight 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


Dimension 


1 





2 





3 


1 


4 


2 


7 


3 


10 


7 


14 


11 


22 


17 



3.2. Approximate conserved quantities. Here we discuss the conserved quan- 
tities of the perturbed KdV equation u% + K(u; e) = 0(e N+1 ) in (2.1). We look for 
the conserved quantity in a formal power series, 



h[u-e] = I Pl (u;e) dx = ^[u] + el^'\u] + ■ ■ ■ + e"l^ ! (u) + 0(e" +1 ), 
so that the density pi(u; e) satisfies 

(3.3) V K ■ pi(u; e) = 0(e N+1 ) mod Imp). 

Here the higher order corrections are given by 



(!)r 



jw[«]= / P f(u)d x , P f> eF 2(l+i)+2 [ 



(i) 



R 



Expanding ( |3.3D in the power of e leads to the equation for p\ 



(m) 



(3.4) 



Then we have, 



Lemma 3.1. For the linear map V. 



V K (0) 



the kernel ojV„(o) with a fixed weight is a one dimensional subspace ofJ^kM given 
by 



KerV^(o) n^2fc[w] = SpanR^ ^}. 

A proof of this Lemma can be found in p3 . Namely there is only one conserved 
density in the form of differential polynomial for each weight. From this lemma, we 
obtain a sufficient condition for the solvability of p^ on the space of differential 
polynomials ^(m+O+a M > 

(3-5) dim f 2 (m+!)+5M = dim -^(m+O+aM - 1- 
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The equation ( |3.3|) is overdetermined in general, and we denote by the number 
of the constraints for the existence of pf""\ that is, 

N^ m) = dim f 2 (m+!)+5M - ( dim •? 7 2(m+0+2 M _ !)• 

Note here that 7V/" ll) = iv/" l2) if /i + mi = Z 2 + m 2 . Then from the Table 1, we 
obtain 

Proposition 3.1. For the existence of the higher order corrections of the conserved 
quantities, we have 

i) There always exist conserved quantities, Iq, Ii, I2, 13 up to order e. 

ii) At order e 2 , there exist 1$, Ii, I2, and = I, that is, there is one condition 
4 — for the existence of I3 where 

^ : = -560af } + 170a ( 2 2) - 60a^ 2) - 8af } + 24a^ 2) - 9a^ 2) + 3a^ 2) + 24a^ 2) 
+f aW(240a« - 674 1} + 6a£ x) ) + |4 1) (4aW + 3o4 1} + a«) 
-4 1) (24 1) +a|), 

iii) order e 3 ? i/iere always exist Iq^Ix. There are total three conditions, ji^ — 

(3) (3) 

0, A; = 1,2,3, luztfi = for the existence of I2 (i.e. iV^ = 1), and two 
conditions fi 2 ' = 0, £i 3 — /or i/te existence of I3 (i.e. — 2). (Explicit 

form of Hk , & — 1)2,3, are listed in Appendix.) 

iv) >lt order e 4 , there always exist Iq. There are total seven conditions, 4 ~ 
k = 1, • • • ,7 mf/t 7V 1 (4) + iV^ 4) + AT 3 (4) = 1 + 2 + 4 = 7. 

One should note here that many physical examples have several conserved quan- 
tities, such as a total mass, momentum and energy, which may be assigned as the 
first three quantities ij[u;e] with I = 0, 1 and 2. Then the existence of the higher 
conservated quantity ^3 [it; e] is a key for the integrability of the perturbed equa- 



tion. The item i) in Proposition 3.1 suggests the (asymptotic) integrability of the 

perturbed equation ( |2.ip up to order e. In fact, in Section g we transform the 

perturbed equation to an integrable system up to order e, and discuss the effect of 
(2) (2) 

the nonexistence of 1^ , i.e. fi\ 7^ 0, on the interaction behavior of two solitons. 

3.3. TV soliton solution. One of the most important aspects of the existence of 
conserved quantities in the form of differential polynomial is to provide several 
exact solutions of the system. For example, ./V-soliton solution can be obtained by 
the variational equation (see for example J24|), 

(3.6) V6{Ai,--- ,Ajv} = 0, 

where 6{Ai, ■ • ■ , An} is the invariants given by 

6 = l£li M + XJ^ [«]+••• + X N l[ 0) [u] . 

Here Afc are real constants (the Lagrange multipliers) . Thus the TV-soliton solution 
is given by a stationary point of the surface defined by [it] =constant subject 

to the constraints, ijP' [it] =constant for k = 1, 2, • • • , N. This characterization of 
the iV-soliton solution is essencial for its Lyapunov stability. 
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The equation (3.6) gives a 2N order differential equation for u{x, ■), and contains 



one soliton solution in the form, 

(3.7) u(x, •) = 2«; 2 sech 2 (K(a; — xo(-)))> 

with a appropriate choice of the parameters In fact, using the bi-hamiltonian 



relation (3.2) the variational equation V6 = can be written as 

(EH^ 1 + mD\ N - 2 + ■■■ + ^_!)D(V4 0) + MwV/f } ) = 0, 
where ^fc's are some constants related to \j. Then by c hoosing //jv = —4k 2 this 



equation admits the one soliton solution in the form (3/7). Here the operator ER is 
called the recursion operator defined as 

(3.8) m=QD~ 1 , and SRDVif 5 = DVlj® v 

The recursion operator plays an important role for the theory of the integrable 
systems. 

4. Symmetry and the perturbed equation 

Here we review the symmetries of the KdV equation and discuss how these 
symmetries work for the analysis of the nearly integrable systems. 

4.1. Symmetries of the KdV equation. Let us define the symmetry for a sys- 
tem, 

(4.1) u t = K(u), with K(u) G V[u}. 



Definition 4.1. A function S(u) G V[u] is a symmetry of (4.1) if S(u) satisfies the 
commutation relation 

ad K ■ S(u) := [K, S]{u) = (V S -K-V K - S)(u) = 0, 

where ad/<- : V[u] — > P[u] is the adjoint map of K . 

This means that if S(u) is a symmetry of K{u) then the flows of vector fields 
generated by K(u) and S(u) commute with each other, i.e. Vk ° Vs — Vs ° Vk = 0. 

The KdV equation has an infinite number of commuting symmetries which are 
given by the hamiltonian flows generated by the conserved quantities J, [u] for 
I = 1, 2, • • • , that is, 

K< n \u)~DVA, for n= -1,0,1,2,- •• . 



Then using the recursion operator EH in (3.8) those symmetries can be constructed 
as, 

K^- 1 \u)=9i n K < f 1) {u), with K ( - 1] {u) = u x . 

The symmetry can be also constructed by the so-called master symmetry |Q. 
The definition of the master symmetries is as follows. A function M(u) is a master 



symmetry of the evolution equation (4.1) if the Lie bracket defined by M maps 
symmetries onto symmetries, i.e., 

[M, Si] = Sj, 
where Si and Sj are the symmetries. 



10 



HIRAOKA AND KODAMA 



The explicit form of the first few commutative symmetries , i > — 1, are 

K { - 1] - u 

1V Q — a Xi 

Kq — 6uu x + u 3x , the KdV equation 
K^ 1 ' = u 5x + I0u 3x u + 2Qu 2x u x + 30u x u 2 , 

K^ = u 7x + I4u 5x u + A2u 4x u x + 70u 3x u 2 + 7Qu 3x u 2x + 280u 2x u x u 
+U0u x u 3 + 70u 3 x , 

K^ = ug x + I8u7 x u + 72u ex u x + 168u 5x u 2x + 126u 5x u 2 + 252u ix u 3x 

,2 i /ion,,„ o,3 



+756u4 X u x u + 1260u 3x u 2x u + 966u 3x u x + 420u 3x u 
hx 1 



-1302uLu x + 2520u 2x u x u 2 + 1260m^u + 630u x u A , 



and the master symmetries Mq , j > 0, are 
M^ 0) = xK^ X) + 2u, 

M (1) = xK { Q 0) + 8u 2 + Au 2x + 2K < f 1) D- 1 {u), 
M (2) = xK^ +32u 3 + A8uu 2x + 36u 2 x + 6u ix + 2K i ° ) D- 1 (u) 
+6K i - 1) D- 1 (u 2 ), 

Af (3) = xK (2) + 8u 6x + 96uu ix + 2A0u x u 3x + 160u^ + 38Au 2 u 2x + 576uu 2 x 
+ I28u 4 + 2K i 1) D- 1 (u) + GK^D-^u 2 ) + 10K { ~ 1] ' D" 1 (2u 3 - u 2 x ). 

The set of those polynomials {K q U \ M^ 1 } forms an infinite dimensional Lie 
algebra containing a classical Virasoro algebra of the master symmetries, that is, 

(4.2) [M n , K m ] = (2m + 3)K n+m , [M n , M m ] = 2(m - n)M n+m . 

One note here that the symmetries Kq (it) are the odd weight (differential) polyno- 
mials, while the master symmetries are the even weight polynomials. In particular, 
each M^ k \u) ~ xK ( k ~ 1] (u) can be considered as an element of the space of even 
weight polynomials generated by the recursion formula starting from k = 0, 



y(*)(u) = a [ k) 



U 2 kx 



+ J2 (4 fe) ( y / fel)Y f 2) )H+^ fe) ( x f 1_1)£, " ly f 2) )( u ))> 



k i +fc2 — — 1 
1 < i < N ( k l) 
l<j<N(k 2 ) 

where X^ k \u) is a monomial in T^fc+sjw]. This defines the space PevenMj 

oo 

^cvcn [u] = V { c k v L [«] , with [u] = pW, [u] ® Q ^ n [«] , 

fe=0 

where 7-evenM = T > 2 k+ 2 \ii\. The space 'PevinM will be important for the normal 
form theory discussed in the next section. 

Remark 4.2. The actions of the symmetries and the master symmetries on one- 
soliton solution (3.7) give a representation of the algebra (4.2) in terms of the 
soliton paramet 
the vector field, 



(n) 

soliton parameters K and 9 = kxq: For the symmetry Kq , the action generates 



V n = 4 n+1 K 2n+3 for n>-l, 
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and for the master symmetry, 



W m =4 n K 2n+1 ^-, for m>0. 

OK 



Those gives a representation of the algebra (4.2), that is, 

Wn,V m ] = 0, [W n ,V m ] = (2m + 3)V n+m , [W n ,W m ] = 2(n - m)W n+m . 
4.2. Approximate symmetries. As we discussed on the approximate conserved 



quantity for the perturbed equation (2.1) in the previous section, one can also 
discuss the approximate symmetries for ( |2.1[ ) . We say that a differential polynomial 



H(u) £ V[u] is an approximate symmetry of order n of the perturbed equation (2.1) 
if the commutator of K(u) and H (u) gives 

[K,H]{u) = 0{e n+1 ). 

Expanding H(u) in the power series of e, H — + eH^ + e 2 H^ + • • • , we have 
the equation for H^ m ^(u), 

m 

ad„ ( o> -if (m) = -y^[KV),H( m -% for m = l,2,---, 

where = u^ x + 6uu x , and K^'a are the higher order corrections of the KdV 
equation ( |2.l| ). Then choosing be one of the symmetries of the KdV equation, 
say H(°) = Kq U \ we find the obstacles for the existence of higher order corrections 
}{i m ) f or the approximate symmetry, and obtain the same conditions for the exis- 



tence as stated in Proposition p.l| |2l| . In the proof, one needs the following Lemma 

1 v n 



similar to Lemma 3.1 for the kernel of the adjoint action ad^to) . 
Lemma 4.1. The kernel o/ad„(o> on the space of differential polynomials P 2n +5H 

Kg 

is given by 

Ker (ad x (o)) n Pan+sM = Span R {i^ n )}, for n > -1. 

5. Normal form theory 

The normal form theory has been well developed in the study of finite dimen- 
sional dynamical systems (cf.Q). The main purpose of the normal form is to classify 
the vector fields near critical points in terms of the symmetries of the leading order 
equation. This concept has been applied for the perturbed KdV equation as well 
as the nonlinear Schrodinger equation in (l4[ [l5| El], [l8) . 

In this section we review the normal form theory for the perturbed KdV equation 
|l5| |l6[| , and discuss the effect of the obstacles on two-soliton interactions. We also 
discuss the Gardner-Miura transformation (which is an invertible version of the 
Miura transformation) in terms of the normal form theory, and show that the 
Gardner-Miura transformation is just a normal form transformation. 

5.1. Normal form. The basic idea of the normal form for the perturbed KdV 
equation is to remove all the nonresonant (nonsecular) terms in the higher order 
corrections using a near identity transformation given by the Lie transformation. 
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Since the symmetries -Kg (u) of the KdV equation give the obvious resonant terms, 
we write each higher order term K^ n '(u) in the form, 

K^\u) = a^ ) K^\u) + F^ n \u), for n = l,2,---, 

so that (u) has no linear term. Then the point of the transformation is to sim- 
plify the term by removing the nonresonant terms in If one succeeded to 
remove the entire F^ n \u) up to n — N, then the perturbed equation is asymptoti- 
cally integrable up to the order e N , and it possesses approximate integrals 7j [v; e] 



for all I £ Z>o and n = 1, ■ • ■ ,N. However as we have shown in Proposition 3.1 



there are conditions for the existence of approximate conserved quantities. Thus 
we expect to see obstacles in removing all the nonlinear terms F^ n \u) from the 
higher order corrections. In fact we have 

Lemma 5.1. There exists a near identity transformation, T e : v i— » u, 
u = T t (v) = u + e0« (u) + ■ ■ • , with ^ (u) S V^Ju] 



such that the perturbed equation (2.1) is transformed to 

(5.1) v t + G{v;e) = 0{e N+l ), 

with G(u;e) = K^°\v) + eG {1 \v) + e 2 G (2) (t>) + ■■■ + e N G m (v), 

where G^ n \u) are given by 

A(n) 

GW(») = a^K^iv) + R [n) {v), with R^(u) = £ ^R\ n \u). 

8=1 

(n) I — ~X 

Here the constants n\ are given in Proposition 3.1 for the existence of the ap- 
proximate conserved quantities, A(n) is the total number of the conditions for the 
existence, and some R^(u) G T^n+sM- 

Proof. We take the Lie (exponential) transform for a near identity transform, 
that is, 

u = T e {v) = cxpV<f, ■ v, 
where the generating function </> is expanded in the power series of e, 
<j> = 60« + e 2 0( 2) + • • • + e W (W) + 0(e N+1 ). 



Substituting this into (2.1) leads to 

Vg = Ad CX p V</) • Vk ■= expl^ ■ V K ■ exp(-V^), 

which gives 

oo 

G = E ^ ( ac W" -K = K+[^K] + -& [0, K}\ + - [</>, [0, [0, K]]] + 

n=0 
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Then from each order of e, we obtain 

[#<j°>,0< 2 )] = kV>-GW + ±[4>W,KW + GW], 



K<$-G<® + $[<f>W,KW + GW], 



Thus we have the equation for <p n > in the form called the homological equation, 
(5.2) ad„ ( o) -ct> {n) := [i^ 0) , = -G (n) , 



where K^ n ' is successively determined from the previous equations. Since the ad- 
action ad„(o) raises the weight by three, we have 



ad„ (0 ) : V 2n +2[v} — > V-2 n +h[v]. 

Then we assume 

0W = a{ 1) u 2 + a { 2 1) u 2x + a i 3 1) u x D- 1 {u), 

J.(2*\ (2) S , (2) (2) 2 . (2) 

<T > = al'u^ + a 2 'uv,2x + a\'u x + a\'v,Ax 



+a { pK$ ) D-\u) + o% ) u x D- l {u 2 ), 



0( 3 ) 



„( 3 )„,4x 



..(3) 



„(3)„,2 



(3)„,2 



„(3)„ 



a\ u + a 2 uu x + ag u u 2x + a\ u 2x + a 5 u x uz x 



+aQUU4 X + a ( 7 'u 6x + a ( 8 'u x D ^u 3 ) + a ( g ' u x D 1 (u 2 x ) 

e 10 ^0 



+a^K^D- l {u 2 ) + a^K^D-^u) + a { ? 2 ] D~ l {u) . 



The homological equation (p.2|) gives a linear system of equations for the column 



vector a 



(5.3) 



in) 



I («) 



' a jV(n)) T Witn = dim ^2«+2M, 



A, 



where A is a (M(n) — 1) x N(n) matrix representation of the linear map ad 

on V2n+2[v], and &(") represents the coefficients of the monomial in — 
which has M(n) — 1 elements with M(n) = dim T^n+sM- The system ( jslil ) is 
overdetermined, and the total number of constraints for the consistency of the 
system is given by 

N(n) := dim V 2n +5 - 1 - dim V 2n +2- 
This number should agree with that of the conditions for the existence of approxi- 



(n) 

mate symmetries, that is, the number of n\ in Proposition 3.1 

2: 



Let us give an explicit form of G ^ n > up to n 
At order e, the matrix A in (5.3) is given by a 3 x 3 matrix with rank 3, so that 
we have R^ 1 ' — 0, that is, no obstacle. The explicit transformation 0( 1 )(u) is given 

by 



„(i) 



(20a\ 



(i) 



,(i) 



4 X) )- 



J 1 ) 



12 
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(i) 



,(i) 



,(ih 



14 



HIRAOKA AND KODAMA 



At order e 2 , the matrix A is an 7 x 6 matrix with rank 6, and we have one 
obstacle with the form 

i? (2) (v) = bf\ 5x v + b^v 4x v x + bfK 3x v 2 + b ( ^v 3x v 2x 
+ bfv 2x v x v + b^v x v 3 + bjpv*, 

(2) 

Here b\ , i = 2, 3, • • • ,8, are constants satisfying the condition, 

170^ 2) - 604 2 } - 8&i 2) + 2Abf ] - 9b { ^ + 36^ 2) + 2Abf ] = 1. 

The explicit formula of the — \ • ■ • , Oi^) T is given in Appendix. □ 

The perturbed equation (2.1) may have several approximate conserved quantities 
based on the original physical setting. Then we consider a particular form of the 
transformed equation ( |5.l|) whose conserved quantities are given by those of the 
KdV equation. We call this form of equation the normal form of (|2.1|), that is, we 
define 

Definition 5.1. For a subset of integers T C Z>o, suppose that the perturbed 
KdV equation (2.1) has the approximate conserved quantities //[w;e] for / G T up 
to order e . Then the normal form of (2.1) is defined by (5.1), 

N 

v t + J2^ n G {n) (v) = 0{e N+1 ), 
n=0 

whose conserved quantities Ji[v; e] := Ii[T e (v); e] for I € T are expressed in terms of 
the conserved quantities of the KdV equation 1^ [v] , 

(5.4) J,[«; e] = 4%] + ec«^M + • • • + e N c^C + M + 0(e N ^), 

where cf\ i = 1, 2, • • • , N, are some real constants. 

In particular, if the set T contains the first three numbers, T D {0, 1,2}, then 
the normal form admits a solitary wave solution in the form of KdV soliton (3.7). 
This can be seen by taking the variation, 

V(J 2 [«;e] + AJi[v;e]) = 0. 

Since many physically interesting systems possess those conserved quantities as 
mass, momentum and energy, we expect to find a solitary wave close to the KdV 
soliton for such systems. Now we show the existence of such normal form for the 
case T = {0,1,2}: 



Proposition 5.1. Suppose that the perturbed KdV equation (2.1) has the first three 
approximate conserved quantities, Ii[u;e], I = 0,1,2 up to order e 3 . Then the 
corresponding normal form takes the form (5.1) with 

RW(v) = 0, 

where the conserved quantities Ji[v,e] for the normal form are given by 
Mv-e] = h[v] + ecfhi+M + ■■■ + e 3 c[ 3) / ;+3 H + 0(e 4 ). 
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The obstacles JZt are expressed as 
Rp = jfc (v 5x v + \v± x v x + bv 3x v 2 - \v 3x v 2x + 20v 2x v x v + 10v x v 3 + 5v 3 ) , 

Ri = (v7 X v + ^v 6x v x + Uv 5x v 2 - \v ix v 3x + 63v 4x v x v + 35v 3x v 2x v 
+56v 3x vl + A2v 3x v 3 + 42v% x v x + 273v 2x v x v 2 + I68v x v 
+ ±fv x v 4 + 2lv x D- 1 (v 3 x )), 

i?2 3) = 1^5 (^xV + I0v 6x v x + ^-v 5x v 2 - I4v ix v 3x + ^-v ix v x v 

+ ^v 3x v 2x v + 203v 3x vl + U7v 3x v 3 + ^-v\ x v x + 903v 2x v x v 2 
+483v x v + 210v x v 4 + 2\v x D~ x (v 3 )) . 

and 9\ is the recursion operator. 

Proof. Recall that J;[v,e] = J R p(v;e)dx is an approximate conserved quantity 
for Q if 

V G -p(v;e) + 0(e N+1 )elm(D). 
Then using the form (5.1) and pi = p[ ' + eCj pf, \ + • • • , we have, at order e 2 , 
(5.5) V R{2) ■ p\°\v) G Imp), for I = 0, 1, 2, 

and at order e 3 



(5.6) 



Vrw ■ Pk\v) 6 Imp), for k = 0, 1, 
V a w -P^+c^V^ -pf Gimp)- 



Then from a direct computation with the explicit form of G 'Pg[w], the condi- 
tions (5.5) lead to the required form of R( 2 \v). 

For (|5.6|), we first write i?' 3 ' in the sum of homogeneous solution R^ and 



a 



(3) 

particular solution R p , that is, 



V r(3) • p[ 0) G Imp), for Z = 0,1,2 
V R l) -fP+cPVum -pf Gimp). 

Then one can find a particular solution by using the recursion operator £H = 0D -1 , 
Vr(2)-4 0) = fl^V/f 3 mod Imp) 

= rWd- 1 dv4 0) mod Im P) 

= R^D^QVI^ mod Im(£>) 
= (GD-^^^V/f 5 mod Imp) 
= V^^) ■ m °d Imp) 



from which wc have 



if) = -p^c^m(R^). 



Then from a direct computation we have the homogeneous solution in the desired 
form, 

4 3) = /4 3) iif +4 3) 4 3) - 
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Remark 5.2. We have the following remarks on the obstacles: 

a) In general, the explicit form of the obstacles W- n > may be successively obtained 
by solving the linear equations, 

n— 2 

V R M ■ P [ 0) + ]T c[ k) V Rl n-v p\% £ lm(D), for n = 2, 3, • ■ • . 

k=l 

Suppose we found up to k = n — 1. Then we have a particular solution 
in the form, 

n-2 

fe=i 

Here one has to check the compatibility among the different I = 0,1, 2. 

b) Since the normal form with T — {0, 1, 2} admits the solitary wave in the form 
of the KdV soliton, one can see that all the obstacles vanish when v assumes 
the KdV soliton solution, 

R (n \v) = 0, when v = 2 K 2 sech 2 (k(x - x )). 

5.2. The Gardner-Miura transformation. It was found in |22| that there is an 
invertible transformation between the KdV equation and the KdV equation with a 
cubic nonlinear term, 

(5.7) u t + 6uu x + — eau 2 u x , 

where a is an arbitrary constant. The transformation is called the Gardner-Miura 
transformation which is an invertible version of the Miura transformation, 

(5.8) v — u — ae x ^ 2 u x — a 2 eu 2 , with a = — \/ a/6. 



Here we treat (5.7) as an example of the perturbed KdV equation, and give an 
explicit formulation of the Gardner-Miura transformation in terms of the normal 
form theory. Since the perturbed equation ( |5.7| ) has an infinite number of conserved 
densities and there is no resonant term as the symmetry of the KdV equation in 
the higher order, the normal form is just the KdV equation. Then we construct the 
normal form transformation, u — T e (v) which is the inverse of the Gardner-Miura 
transformation (5. J). 

Since the Gardner-Miura transformation (5.8) is of a Riccati type, the change of 
the variable 

(5.9) u= - ( D]n<p+ —) , with 5 = ae 1 / 2 , 

o \ 25 J 

leads to the Schrodinger equation, 

L 2 <p := (D 2 + v)(p = k 2 cp, with k= -1/(25). 

Using the notion of the pseudo-differential operators D v for v G Z, one can define 
L = (D 2 + v) 1/2 as 

L = D + c h D- 1 +q 2 D- 2 + ••• , qi eV i+ 2[v], 

Then writing D in the power series of L, we have 

D = L+p 1 L- 1 +p 2 L- 2 + -- - , Pl eP l+2 [v] 
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and using Lip — kip, we find 

^ oo 

(5.10) DU<p = -—+Y f (-28) i p i (v). 



28 

Note in particular that we have 



From ( |5.9| ), the equation ( 5.10 ) leads to the inverse of the Gardner-Miura transfor- 
mation, 

oo 

u = ^(-2)M i -Vi(v,Vx,---) 
»=1 

= w + 8v x + 6 2 (v 2x — v 2 ) H . 

We now remove the terms of the non-integer powers of e, the odd integer of 8 = 
ae 1 / 2 , in this equation. The first term 8v x can be removed by the translation of 
x. After removing the term 8v x , we have Kq(v) at the order of e 3 / 2 . This can be 
removed by the translation of t. Continuing this process, we see the symmetries 
of the KdV equation at the non-integer powers of e. Then shifting the symmetry 
parameters £2,1+1, those can be removed. Here x = X\ and t = —X3. Now we can 
show 

Proposition 5.2. The inverse of the Gardner-Miura transformation gives a nor- 
mal form transform, 

u — -^D^ 1 sinh ( — — - — ] v, 8=—\Mr, 

S \ ^ 2i + 1 dx 2t+1 J V 6 

where the derivative of v with respect to X2n+i defines the symmetry, that is, 
dv/dx 2n +i = K { n ~ 1] {v). 

Proof. It is well known |9| that the wave function ip of the Schrodinger equation 
can be expressed by the r-function, 

^(x )fe ) = li^Me^, with fc = -i, 
t(x) 2d 

where we denote x := {x\, x 2 , ■ ■ ■ ), and 

S 3 6 5 



W = (S'5 

Expanding the equation D In ip with this equation in the power of 5, we find 
(5.11) u(x) = 2L> 2 lnr(x). 



Then from (5.9) we have 
Now applying the vertex operator 



u(x) = - C D (lnr(x + 2(8}) - lnr(x)) . 




I 00 ij2i 



(2i + 1) dx 2 i+i 
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and using (5.11), we obtain the result. Note here that both u(x) and u(x + (5)) 
satisfy the same equation (5.7). □ 

6. Interactions of solitary waves 

As an important application of the normal form theory, we consider the interac- 
tion problem of solitary waves, and show how the theory enables us to understand 
the interaction prop erties under the influence of the higher orders in the perturbed 
KdV equation ( |2.l|) . We assume that the perturbed equation possesses the first 
three conserved quantities, i.e. T — {0, 1, 2}. Then the first obstacle appears in the 
order e 2 as fj,^ ^ 0, that is, we consider the normal form, 

v t +K^°\v) +ea[ 1) K i 1 \v) + e 2 {a {2) K {2) {v) +^R[ 2 \v)) = 0. 

Several physical examples of this type will be discussed in Section ^. We use the 
method of perturbed inverse scattering transform [ fp2[ [l3| to analyze the normal 
form. We start with a brief description of the method. 

6.1. Inverse scattering transform. The key of the inverse scattering transform 
is based on the one-to-one correspondence for each t between the scattering data 
S(t) and the potential v(x,t), the solution of the normal form, of the Schrodinger 
equation (see for example 0), 

+ (« + * s V = o. 

The correspondence is given by the formula of v(x, t) in terms of the squared eigen- 
functions, 

(6.1) v(x,t) = 4 Kj(t)Cj{t)<p 2 (x,t;iKj) + — / kr(t;k)<p 2 (x,t;k) dk. 

Here the scattering data S(t) is defined by 

S(t) := {{Kj(t) > O.CfC*)}^!, r{t]k) for k E R} , 

and the eigenfunction ip(x, t; k) is assumed to satisfy the boundary condition, 

tp(x, t; k) — > e~ %kx , as x — > — cxo. 

Note that the squared function ip 2 (x,t; k) satisfies 

Qip 2 = (D 3 + 2(Dv + vD))ip 2 = -4k 2 D(p 2 , 

so that the function Dip 2 is the eigenfunction of the recursion operator 9t with the 
eigenvalue — 4fc 2 . 

With the formula of v in (3.1), the conserved quantities Ji [v ; e] can be expressed 
in terms of the scattering data, 

(6.2) J t [v; e] = /| 0) [v] + ec^I^v] + ^cfhfhv] + 0(e 3 ), i = 0, 1, 2, 



where /• can be expressed as 



o2m+l 

(6.3) 4°) = ^-—^-ir^ 



M 

X>? m+1 - (-l) m A, 
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with the radiation part A m (t), 

A m = -^+^ [°° k 2m \n(l-\r(t;k)\ 2 ) dk > 0. 
2tt Jo 

The existence of such conserved quantities plays a crucial rule for the interaction 
mechanism as in the case of the KdV solitons. 

The time evolution of the scattering data is determined by 

£-S(t) = -V K > S(t), for v t + K(v) = 0. 
dt 

In particular, the equation for the reflection coefficient r(t; k) is given by 
d 1 f 

dt r[t]k) = 2ika>(t;k) /^fr*)^.*;*) dx. 

where a(t; k) is the reciprocal transmission coefficient determined by r(t; k) and Kj. 
Using the normal form vt + K(v) = 0, we obtain the equation for r(t; k), 



dr _ . 2 $ } 



R (2) ip 2 (x,k) dx + 0(e 3 ), 



dt 2ika 2 

where to = 8fc 3 (l - iea^k 2 + 16a^ 2) e 2 fc 4 ). 

Those given above provide an enough information for our purpose of studying 
the interaction of solitary waves. 

6.2. solitary wave interaction. Recall that the obstacle vanishes for one soliton 
solution of the KdV equation, that is, there is no effect of the obstacle on the 
solitary wave. The higher order terms lead to the shift of the velocity of the soliton 
solution due to the resonance caused by the symmetries of the KdV equation jl9| . 
In order to see the effect of obstacles, we now consider the interaction of two solitary 
waves, and show the inelasticity in the interaction which can be considered as a 
nonintegrable effect of the obstacle. 

6.2.1. Inelasticity in the interaction. Let us assume that the initial data consists 



of two well-separated solitary waves with parameters K\ > K2 in the form of (3.7) 



traveling with speed Sj = 4/t5-(l — Aea\ L ' k 2 + I6e 2 a^ k*) j — 1, 2, and approaching 
each other. Then we analyze the interaction by using a perturbation method where 
the leading order solution is assumed to be the exact two-soliton solution of the 
KdV equation, that is, for large xqi — X02 ^> 1 

v(x, 0) ~ 2K^sech 2 (Ki(a: — Xoi)) + 2K 2 ,sech 2 (K2(2; ~ ^02))- 
We then determine the evolution of the scattering data, Kj (t) and the radiations 



A m (i) in (3.4). First we have 



Proposition 6.1. |lq ] Due to the interaction of two solitary waves, their param- 
eters Kj j — 1,2 are shifted by AKj(t) which can be expressed in terms of the 
radiations A m (t), 

5k|Ai+3A 2 n , A . n . 5k?Ai+3A 2 ^, A , „ 

Aki - iskk^i + 0(eAm) - °' AK2 = - 15^-^) + ° {eAm) - °- 

There is also a production of the third soliton with the parameter, 

A« 3 = 15 ^ A ° + 5 H + 2 Ki)Al+3A2 + 0(eA m ) > 0. 
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Proof. From the conserved quantities Ji [v; e] ^constant for I — 0, 1, 2 up to order 



e 2 , we obtain, 



AI, 



(0) 



0(e 3 



where the variations AJ, 
that is, 



(0) 



are taken over the shifts Ak, and the radiation A, 



(21 



-l)'Ai = 0(eA n ), for 2 = 0,1,2. 



JV 

l)5>?A/e 3 -- 
i=i 

Here iV — 2 is a possible number of new solitary waves. Since the Anj for j ' > 2 
represents a new eigenvalue of the Schrodinger equation, it is nondegenerate and an 
isolated point on the imaginary axis of the spectral domain . This implies that the 
number of new eigenvalues should be just one for a sufficiently small e, i.e. N = 3. 

Because of Kj = initially for j > 2, one can find the formulae of Akj for j = 1,2 
from those variations for I = 1, 2. Also the first variation with 2 = 0, 

3 

^A«y-A = O(eA„), 
i=i 

leads to the formula A«3 for a new solitary wave. □ 
Thus we find the followings: 

a) The total mass, M = J vdx oc k, of the larger solitary wave is increased, and 
contrary that of the smaller solitary wave is decreased by the interaction. 

b) The amount of the energy change, E = J v 2 dx oc k 3 , has the property 1 < 
\AE 2 /AE-l\ < (k 1 /k 2 ) 2 , i.e. the energy expense of the smaller solitary wave 
is more than the energy gain of the larger one. 

c) The interaction produces a new solitary wave as well as radiation. 

Those results except for a new solitary wave production are consistent with the 
numerical observation in Q. It may be difficult to observe the new solitary wave 
from the numerical calculation, since this solitary wave has long width and small 
amplitude (Ak^) 2 which is of order e 8 (see below). 

The function A m (t) of the radiation can be computed as follows: The reflection 
coefficient r(t; k) in ( |6.l[ ) can be expressed as 



r(f ; k) 



2ik 

Then A m (t) is given by 
(6.4) A m (t) 
where D m (t; k) is 



D m (t;k) = k 2m \r{t;k)\ 2 



dr- 



a(r; fc) 5 



dxLp 2 (x,T]k)R { ?\v) 



2m + 1 



: (/4 2) ) 2 fc 2(m ~ 



2tt 



o(e 4 ) 



D m {t;k) dk + o(e 4 ), 



dr e 



dx ip 2 (x,T;k)R^(v) 



o(e 4 ). 



We numerically calculate the formula A m (t) by means of perturbation, that is, we 
assume v(x, t) to be a two-soliton solution of the integrable one v t +ea^ Ki + 

and (p(x, t, k) is the corresponding eigenfunction. The result will be 



e 2 a\ 2) K 2 
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shown in Section 7 for the examples of ion acoustic waves and the Boussinesq 
equation. 



6.2.2. The additional phase shifts of solitary wave. As a consequence of the nonlocal 
terms in the normal form transformation, one can find the additional phase shifts 



on the solitary waves u{x, t) of the perturbed KdV equation (2.1) through their 
interaction. 

Let us first recall the phase shifts of the two-soliton solution for v(x, t) ]8|, g|. 
The asymptotic form of v{x, t) consists of well separated one solitons, 

± f t —>■ ±oo 

v{x,t) « u, (x,t), as < 

J [ KjX ~ Wjl 

with v^(x, t) — 2K 2 sec\v{KjX — uijt — Oj). 
Then the phase shifts Axj := [fff — Oj)/nj are given by 

Ax ( V = -Un( Kl - K \ ) and Ax^ = -1 In f Kl = " 2 

We now compute the correction to the shift Ax^ using the normal form trans- 
formation. From the asymptotic form of the two-soliton solution for t — > — oo, we 
have up to order e 

+ea 3 (?; 2 ") a; / dx. 

One should note here that there is an extra term in the solitary wave which is 
the key term for the additional phase shift. Namely the term can be absorbed as 
a translation of x in vi. The other terms contribute to modify the shape of the 
soliton, the dressing part. 
Also for t — > +oo, we have 

u(x,t) — » v+ + e(a[ 1) (v+) 2 +a { 2 1) (v+) 2x +a^(v+) x f v+ dx 

rei x~wi« \ - J_ OQ 

+ea 3 (v^) x J dx. 



u(x,t) — > v 1 +e[a\'(v 1 ) + a\ (ui )ix + <A (% ) x / v 1 dx 

K\X^UJ\t 1 ' 



— ► D 2 +e a '("2) + a 2 ( v 2 hx +ay(v 2 ) x / w 2 dx 



e ( a« («+) 2 + a« (« 2 + )2, + 4 X) («2 + ). / «2 + <** 



Now the additional shift appears to v\. This can be extended for the next order 
where the shift also appears as a traslation of t with a term like Kq°\vi) J R V2dx. 
Thus we have the total phase shift for v\ solitary wave, 

(6.5) Arcx = Axf ] + eAx[ 1} + e 2 Ax[ 2) + 0(e 3 ) 
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where the additional phase shifts are given by 



Axi = / V2<lx = An^a^p , 



7. Examples 

In this section, the normal form theory is applied to some explicit models in- 
cluding the ion acoustic wave equation, the Bousinesq equation as a model of the 
shallow water waves, and the regularized long wave equation (sometimes called 
BBM equation). We also carry out the numerical simulation for those examples 
and compare the results with the predictions obtained from the normal form theory 



such as the phase shift ( p.q ) and the radiation energy (6.4). 

We also consider the 7th order Hirota KdV equation which is known to be 
nonintegrable even though it admits an exact two solitary wave solution. The main 
issue is to determine the order of the obstacle of the corresponding normal form, 
which indicates an nonintegrability of the equation in the asymptotic sense. It 
turns out that the obstacles appear at order e . 



7.1. Ion acoustic waves. An asymptotic property of the ion acoustic waves has 
been discussed in several papers (see for example j2§]). These studies show that 
the KdV equation is derived as the first approximation of the ion acoustic wave 
equation under the weakly dispersive limit. The higher order corrections to the 
KdV soliton solution have been also discussed in [pl| . Recently, Li and Sattinger in 
|20[ studied numerically the interaction problem of solitary waves and showed that 
the amplitude of the radiation after two solitary wave interaction can be observed as 
small as 10 ~ 5 order, and they concluded that the KdV equation gives an excellent 
approximation. Here we explain those observations based on the normal form 
theory developed in the previous sections. 

Following the method in |^8| , we first derive the KdV equation with the higher 
order corrections. The ion acoustic wave equation is expressed by the system of 
three partial differential equations in the dimesionless form, 




where n,-, u,-, and <f> are the normalized variables for ion density, ion velocity and 
electric potential. Electron density n e is related with (f> as n e = exp(</>). Assuming 
the weak nonlinearity and the weak dispersion, we introduce the scaled variables, 



vi = ev, n e = 1 + en, and x = e 1 ^ 2 X, t = e 1 ^ 2 T. 
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Then we write (7.1) in the following form for (n,v), 

1 



I -elf 



v t + vv x + 



1 + en 
1 



n t +v x + e(nv) x - e(v ln(l + en)) x = 0, 



0. 



l + en 

Then inverting the operator in front of n t , we have the matrix equation, 



(7.2) 



— U + A o — U + e—B(U) = 0, with U 



1 

1 



where Aq is the constant matrix A® 
given by the expansion 

B = B^ +eB^ +e 2 B^ + ... 



, and the vector function B(U) is 



with < 



BCD = 



B (3) 



v 2x + nv 



Uv 2 ~n 2 ) 



v 4x + n x v a 



5(2) = 



-in 4 



Thus for the case with e = 0, we have two simple linear waves propagating with 
the speeds A± = ±1 given by the eigenvalues of the matrix Aq. We then look for 
an asymptotic wave along with the speed A+ = 1, so that we introduce the scaled 
variables on this moving frame, 

which gives, after dropping the prime on the new variable, 



(7.3) 



dU dB(U) \ 
~dt + dx 



0. 



where I is the 2x2 identity matrix. Let us decompose U in the form, 

U(x,t) = u{x,t)R + + f(x,t)R-, 

where R± are the eigenve ctor s corresponding to the eigenvalues A± = ±1. Then 
taking the projections of (7.3) on the R± -directions, we have 

u t + {L+B) x = 0, 

-2f x + e (/ t + (L-B) x ) = 0, 

where L± are the left eigenvectors with the normalization, L±R± = 1, L±R T = 0. 
Then we can see that f(x, t) can be expressed in an expansion form, 

f(x 1 t)^ef^(u) + e 2 f^(u) + e 3 f^(u) + --- , with /<*>(«) G ?«„[«]. 

where are determined iteratively from 

/ = 6 -(L.B) + e -D-^f t = |(L_B«) + e j(L_B&) + \f t + • • •. 
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Thus we obtain the perturbed KdV equation which takes the following form up to 
order e 2 after an appropriate normalization, 

A 3 ?\ 2 f 5 1 

u t + 6uu x + u 3x + e I -u 5x - u 3x u - -u x u I + e I jjU 7x - -u 5x u 

17 1 3 2 1 7 3 1 3 \ ^. 3 , 

--^U ix U x + -U 3x U 2 x + -^U 3x U - -U 2x U x U - -U x + -U X U I = 0(6 ), 



from which we have 



(2) 

Thus the normal form of the ion acoustic wave equation has the obstacle R\ . The 
normal form transformation is given by 

.A) 



„( 2 ) - 



J 2 ) 



11 (1) 

12' 2 


1 (i) 
3 3 


7 

~ 6' 




731 (2 ) 
3600' 2 


89 
~ 600' 


(2) 

a 3 = 


433 
3600' 


51 (2) 

800' 5 


23 

1800' 


(2) 
"6 = 


87 
400' 



Remark 7.1. The physical variables n and v are expressed as, 

' 1 



n = u + /, w = it — /, with i?± = 



±1 , 



Since / has an expansion of the power series of e and each coefficient f( n > is an 
element in T^evcn; the expressions of ti, v have the same form as of the normal form 
transformation. In the asymptotic sense, all the physical variables are expressed 
by one function u as U = uR + . Then the choice of the higher order terms / has 
a freedom. Then the main purpose of the normal form is to use this freedom to 
classify near integrable systems in the asymptotic sense. 

We now compare the results in Section 6 with numerical results. We used the 
spectral method [ |30"| for the numerical computation. For a convenience, we consider 
( [7.l| ) in a moving frame with a speed c, 

(n) t - (6 + c)n x + 6v x + 6(nv) x = 0, 

= 0, 



(v) t -(6 + c)v x +6\ — +, 

4>2x - 3 exp 4> + 3(n + 1) = 0, 

The computation is done with the 2 12 number of Fourier modes and the time step 
dt = 0.008. In general, the spectral method yields eliasing errors from the high 
frequency modes in the nonlinear terms, so we try to get rid of their errors by 
adopting the 3/2 rule ||. 

Let us first comapare a solitary wave solution of the ion acoustic wave equation 
with the KdV soliton solution. Since the one-soliton solution (3.7) is a kernel of the 
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obstacle R\ , one can construct a solitary wave for the variable u by the normal 
form transformation, i.e., 

u = v + + e 2 U% + \ V, • + 0{e% 

with v — 2n 2 sech 2 {k(x — Xq)), 

where xo is a center position of the soliton at t = 0. The solitary wave for the 
variable u is constructed by the core v with the dressing terms +e 2 (4>^ + ^^< 1 > ■ 
0^^) + O(e 3 ). Now let us see the effect of the dressing terms, which is considered as 
nonresonant higher harmonics in a periodic solution of finite dimensional problem. 
These three figures show the difference of the amplitude of radiations emitted during 



U 

0.002 
0.001 



-0.001 





50 X 30 40 50 X 3° 

Figure 1. Eeffects of the dressing terms. 




50 X 



the time evolution for each initial wave. The initial solitary wave is given by u - 
in the left figure, u = v + e^ 1 ) in the central figure, and u = v + €<j>^ + e 2 (<ft( 2 



') in the right figure. The core part is given by ( |7.1| ) for k — 0.2. These 
figures show that the generation of the radiation can be suppressed by adding the 
higher order dressing terms into the core. 

Now let us discuss the interaction of the two solitary waves. In Figure 2, we 
show the result of the phase shift of the solitary wave with the parameter «i(> K2) 
for various values of K\ — «2- Here we fix k\ + K2 = 0.5 and e = {n\ + k\)/2 ~ 
0.07. The solid line is calculated from the formula ( |6.5[ ), and the broken line is 
obtained by the numerical simulation. The phase shift of the KdV soliton is also 
shown as the dotted line. As we can see that the phase shift formula gives a good 
agreement to the numerical results. The energy of the radiation emitted after the 
interaction is caluculated from (6.4) and is shown in Figure 3. This shows that the 
radiation energy is of order 10 -5 which also agrees with the result in (20). Thus 
the normal form theory provides an accurate description of the deviation from the 
KdV equation. 

7.2. Boussinesq equation. The Boussinesq equation as an approximate equation 
for the shallow water waves is given by 



Vt + v x + (v v )x = 0, 

I, Os 1 



v t + r(« )x + Vx ~ t; v xxt = 0, 
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400 



Figure 2. Phase shift. Figure 3. Time evolution of Ai(t). 



where r/ and v are the normalized variables which represent the amplitude and the 
velocity |^9| . Since this equation is truncated at the first order from the shallow 
water wave equation, the normal form may not provide a structure of the asymp- 
totic integrability at the second order. However the Boussinesq equation itself has 
an interesting mathematical structure, such as the regularization at the higher dis- 
persion regime, and it may be interesting to study its own asymptotic integrability. 
The normal form for the shallow water wave equation has been studied in ]l7j . 

Following the similar process as in the case of ion acoustic waves, we obtain the 
perturbed KdV equation up to order e 2 , 



u t + 6uu x + u 3x + e -u 5x + -u 3x u + Au 2x x - -u x u + 



32 



U7x 



11 
16 



u 5x u 



99 95 3 2 33 27 , 3 , 

-^u ix u x + —u 3x u 2x - -r^u 3x u - —u 2x u x u - —u x + —u x ii 
32 16 16 16 32 16 



= 0(e 3 



("2) 

Due to the non-zero integrability condition ' = 3/2), the obstacle R\ does 
not disappear in the second order correction. The generating functions in the Lie 
transformation are given by 

13 (i) 17 Hi 3 



„(!) 



J 2 * 



a 



24' 



(2) 



2591 
14400' 

1583 



a 



,.(2) 



(2) 



a 

71 

" 100' 

17 

800' 



(i) 



4' 



,(2) 



,.(2) 



2747 
3600' 



571 



9600' 800 4800 
Now let us consider the phase shift during the two solitary wave interactions in the 
same way as the ion acoustic wave equation. The result is shown in Figure 4. 

Again, the numerical data of the phase shift can be well explained by the phase 
shift formula. The time evolution of Ai(t) is shown in Figure 5. This figure shows 
that the energy of the radiation after the interaction is about the same as that for 
the ion acoustic wave equation. 
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Figure 4. Phase shift. 



Figure 5. Time evolution of Ai(i). 



7.3. Regularized long wave equation. Since the KdV equation is derived under 
the assumption of weak dispersion, it is not valid for the wave phenomena involving 
short waves. As a model of shallow water waves for a large range of wavelength, 
Benjamin et al in ||, proposed the equation 

Wt + wx + 6wwx — WXXT = 0. 

This equation is called the regularized long wave equation (RLW) or the BBM equa- 
tion. It has been shown that this equation has only three non-trivial independent 
conserved quantities [p5f indicating its nonintegrability. The BBM equation admits 
a solitary wave solution, 

2k 2 



rsech | kX 



I -An 2 \ I -4k 2 ' 

The numerical simulations by Bona et al. || showed that the interaction of the 
solitary waves is inelastic and generates radiations. Their study also found the 
shifts of the amplitudes of two solitary waves after the collision. The normal form 
theory has been applied to this equation in Jl5| . Here we review the work [ fl5| 
and add the phase shift results. We also compare the BBM equation with the ion 
acoustic wave equation and the Boussinesq equation. 

As in the previous cases, we introduce the scaled variables, 

x = e^(X — T), t = eiT, w = eu. 

which yields 

(1 - eD 2 )u t + 6uu x + u 3x = 0. 
Then inverting the operator in front of Ut, we obtain the perturbed KdV equation, 



OO 

fc=i 



e : D 2k K, 



Then the obstacle R^ 2 \u) appears with fj,^ = 40. The comparison of the phase 
shift between the numerical results and those by the formula (3.5) is shown in 



(2) 

Figure 6. Because of the large value of fi\ , the agreement between the numerical 
results and the results from ( |6.5[ ) is poor. In the computation, the value of is 
about of order one, and thus one needs to consider much higher corrections to get a 
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Figure 6. phase shift. 



better agreement. However the numerical observation of the radiation agrees with 
the formula (3.4) which is of order 1CP 2 ~ 10~ 3 , and the shifts in the parameters 
observed in JBJjalso agrees with the result from the normal form theory. 

7.4. 7th order Hirota KdV equation. Here we consider 

wt + w 7 x + 28ww 5 x + 28wxW4x + 70w 2 xW3x 

(7.4) 

+210w 2 w 3X + A20ww x w 2 x + 420w 3 w x = 0, 

which has a Hirota bilinear form, 

Dx{D 7 x +D t )t-t = 0, 

where Dx and Dt mean the Hirota derivative, and w is given by w = 2D 2 lnr. 
If the order of the derivative is either three or five, the equation become the KdV 
equation or the Sawada-Kotera equation which are both integrable. However ( [7.4[ ) 
is known to be nonintegrable |23| , but it admits both one and two solitary wave 
solutions in the same form as the KdV solitons except their time evolution. Since 

(n) 

there is an exact two-solitary wave solution, several conditions fi k — should 
be satisfied. Now the question is to determine the order in which the condition 

A*!™'' = breaks. 

In order to apply the normal form theory, we first introduce 

w = eu + c, t = 210e 3/2 c 2 T, x = e 1/2 (X - 420c 3 T), 



which puts (7.4) in the form (2.1) of the perturbed KdV equation 

Ut + Kq°\u) + e J^( U Hx + 15 u 3x u + 15u 2x U x + 4:5U X U 2 ) + ^j~2 ( U7x + 2& u 5xU 

+28u4 X u x + 70u 3x u 2x + 210u 3x u 2 + 420u 2x u x u + 420u x u 3 ) = 0, 
where c is an arbitrary non-zero constant. The direct calculation shows 



< 



a4 2) - o, 

l4 a) =0, i= 1,2,3, 
a4V0, z = 1,2,--- ,7. 



Thus the 7th order Hirota KdV equation passes the asymptotic integrability con- 
ditions not only at order e 2 but also at order e 3 , and the first obstacles appear 
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at order e 4 . Then the normal form up to order e of the 7th order Hirota KdV 
equation takes the form, 

(7.5) 

v t + + eaPKP + M 2) 4 2) + ^i } Ko } + e W + i? (4) ) = 0(e% 

where consists of seven obstacles. Since the 7th order Hirota KdV equation 
admits the two-soliton solution, the kernel of R^ should include not only one- 
soliton solution but also two-soliton solution. Thus we have 



Corollary 7.1. If the normal form of a perturbed KdV equation has the form (7.1) 



then there exists two-soliton solution up to order e 4 . Equivalently, if pi{\ /i^ 
solution up to order e 4 . 



(3) (3) 

[i 2 , and /ig are all zero, then the pertubed KdV equation can admit a two-soliton 



Appendix 

Coefficients of the generating function (jP^ : 

(2) _ _loo r (l)s 2 _ 10„(1)_(1) _ _5_( n (V)2 _ ± n W„W + ±(J 1 h? 
a l — 9 \ a l I 27°1 fl 2 108^2 / 54 fl 2 "4 + 36 V a 4 ) 

+ il2 a ( 2 ) , 2 (2) 1 (2) _ 1 (2) _ 4 (2) 
' 9 1 ^ 9 2 ' 9 U 4 6 7 225 ^1 ' 

a 2 ' = — "ir( a i "I" !T a l ' a 2 ' + Ts a l ^ a 3 "* — 3l; a 2 ^ a 3 ' — T8 a l ^ a 4 ' — T8~ a 2 "* a 4 ' 

36 a 3 "4 T 18 V u 4 / ^ 3 "l 2 2 ^ 12 U 6 4 7 75 "l ' 

a ( 2 ) - _2B f «\2 _ JL a (l) a (l) + J_f a «)2 , A a ( 1 ) a ( 1 ) + J- a (1) a (1) - -5-a (1) a (1) 

u, 3 — 3 yu, 1 ) 10 8 u l "2 ^ 108 ^ 2 I T 18 "l u 3 ^ 72 u 2 u 3 36 1 "4 



216 U 2 u 4 24 "3 a 4 ^ 24 V u 4 i + 9 "l ' 18 u 2 18 u 4 6 U 7 

^6 U 8 900^1 ' 

M — 36V U 1 J ' 48 U l 2 144 u 2 u 3 72 u l u 4 288 U 2 u 4 288 3 u 4 

4 .J^C ( 1 )l2 , 7 (2) _ J_ (2) 1 (2) 1 (2) _ j_ (2) _ 7 (2) 
+ 144 l - a 4 7 + 6 U 1 48 U 2 ' 24 u 3 ^ 96 u 6 32 u 7 2400^1 ' 



5 — j l a l J + 9 B 1 tt 2 + 3 "l 3 U 2 150 "l ' 



a 6 — 3 l a l ; 18 a l °2 + 3fV a 2 ; T 18 a l a 4 + 36 a 2 a 4 3 a l 

, 7 (2) _ 1 (2) , 1 (2) 
" l "3 a 2 3 a 4 + 75^*1 • 

The values ^t^, 3 ' for k = 1, 2, 3: 
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2|0 a (3) _ 280 fl (3) + g 0a (3) _ ^(3) _ 5 fl (3) _ fl (3) + ^(3) _ fl (3) _ ^(3) 

+§ag> - ag> - 2ag> + f>(a«)3 + f^) 2 ^ + f^W) 2 - K^) 3 

-f^M 5 - + u^) 2 ^ + fai 1) (4 1) ) 2 + 

— 5?(4 ') 3 + Tf(4 ^4 ^ — 2t4 ^4 '4 ' — 3?4 '4 ^4 ^ ~ Hp4 ^4 

~^TT a 2 ^4 ' ^ fp4 ^4 ' fr4 ^4 ^ + lf4 ^4 "I" IT 4 ^4 ^ + ^ir4 ^4 ^ 

_10-(l)-(2) _ 20-(l)-(2) , 5.(1). (2) _ M^U 2 ) i - IflW,,* 2 ) 

g u 2 u 3 3 u 3 u 3 ^ g a 4 u 3 g a l u 4 ' 3 "2 u 4 3 a 3 u 4 

+a< 1) a< 2) - 12a«4 2) - la^af) + fa^a* 2 ) + 



-42004 3 } + 7704 3) - flaf - 4o4 3) - ifaf + 4o4 3) + f 4 3) - 84 3) 
-6a< 3 > + 7ag> + 8a<»> - ^ + 2ag> + 9ag> - ^o^) 3 + ^R^) 2 ^ 

4V) 2 + #(4 1} ) 3 + ^(4 1} ) 2 4 1} - fa}W 
+1 (4 1) ) 2 4 1) + ?4 1} (4 1} ) 2 - |4 1, (4 1) ) a + £(4 1} ) 3 - f^(4 1) ) 2 4 1) 

+ 10a«4 1) 4 1) + Pa^a^a^ + ^'of - ^a^a™ 

_560„(1)„( 2 ) X 280.(1) -(2) _ 15950 .(1) J2) , 685„(1)„(2) , 310.(1).(2) 
g u 3 f 3 u 4 a >1 g u 2 1" 9 u 2 u 2 ' 9 u 3 u 2 

_6B_(i)_(2) , mJVnW - *°„W„W - ™n W a (2) + 310« (1) fl (2) 
-g-a 4 U 2 i ~ u i "3 IT 2 u 3 IT 4 a 3 ^ ~ 1 4 

-|4 1} 4 2) - 4 1} 4 2) - 2§0 a (D a (2) + 284D42) _ 4^)42) + | a (D a (2) 

+~ir4 ^4 ' — \ a 2 ^4 ' — T a i ^4 ' — ir4 '4 ^ + 14 ^4 ^ 

_4 a (l)J2) _ 7 (2) 
3 a 3 "8 3^1 ' 
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a4 3) = 5604 3) - if^af + 300ai 3) - ipaf - 1004 3) + f af %af + 40<4 3) 
-a< 3 > 20a™ + w (3) - - 2a< 3 4 > - 10a< 3 B > + ^f^) 3 

+i(4 1) ) 2 4 1) + f ai 1) (4 1, ) a - ¥ 4 1} (4 1} ) 2 - ^(4 1} ) 3 - ^(tfW 

+|(4 1) ) 2 4 1) - f 4 1} (4 1} ) 2 - ^a^a^ - 5604 1} 4 2) + ±f>4 1} 4 2) 

+28041)42) _ ifo a w a ( 2 ) + 9o4 1) a ( 2) + f 4 1} 4 2) - ^4 X) 4 2) 

+400ai 1) 4 2) - 200 a (D a (2) _ so a (i) a (2) + 85 a (i) a (2) + M a (D a (2) 
_ 141)^) _ 374D42) + i a (D a (2) _ 160a (D a (2) + 50 a (i) a (2) 

,22 (1) (2) _ 13 (1) (2) _ 85 (1) (2) _ 7 (1) (2) , 5 (1) (2) 
"r" 3 a 3 a 5 3 a 4 a 5 9 a l a 6 2°2 a 6 + 3 a 3 a 6 

— ^~ ^6 ~3~^1 ^7 " ~3^"2 ^7 ^7 _ ^~ ~3~^1 ^8 

+a (1) a (2) + 2 a (1) a (2) - W + 7 -u (2) 

i a 2 a 8 + 3 tt 3 a 8 3 a 4 "8 ' 3 Ml • 
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